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Abstract 



We consider the medium- and long-baseline oscillation physics capabilities of 
intense muon-neutrino and muon-antineutrino beams produced using future 
. upgraded megawatt-scale high-energy proton beams. In particular we con- 

sider the potential of these conventional neutrino "superbeams" for observing 



i-£h ' —* v e oscillations, determining the hierarchy of neutrino mass eigenstates, 

and measuring CP-violation in the lepton sector. The physics capabilities of 
superbeams are explored as a function of the beam energy, baseline, and the 
^Cj . detector parameters (fiducial mass, background rates, and systematic uncer- 

tainties on the backgrounds). The trade-offs between very large detectors with 
poor background rejection and smaller detectors with excellent background 
| rejection are illustrated. We find that, with an aggressive set of detector pa- 

rameters, it may be possible to observe — > v e oscillations with a superbeam 
provided that the amplitude parameter sin 2 26\% is larger than a few xl0~ 3 . 
If sin 2 2#i3 is of order 10 -2 or larger, then the neutrino mass hierarchy can be 
determined in long-baseline experiments, and if in addition the large mixing 
angle MSW solution describes the solar neutrino deficit then there is a small 
region of parameter space within which maximal CP- violation in the lepton 
sector would be observable (with a significance of a few standard deviations) 
in a low-energy medium-baseline experiment. We illustrate our results by 
explicitly considering massive water Cherenkov and liquid argon detectors at 
superbeams with neutrino energies ranging from 1 GeV to 15 GeV, and base- 
lines ranging from 295 km to 9300 km. Finally, we compare the oscillation 
physics prospects at superbeams with the corresponding prospects at neutrino 
factories. The sensitivity at a neutrino factory to CP violation and the neu- 
trino mass hierarchy extends to values of the amplitude parameter sin 2 26\z 
that are one to two orders of magnitude lower than at a superbeam. 



I. INTRODUCTION 



Measurements MM of the neutrino flux produced by cosmic ray interactions in the at- 
mosphere || have led to a major breakthrough in our understanding of the fundamental 
properties of neutrinos. The early observations of the atmospheric neutrino interaction 
rates found a muon-to-electron event ratio of about 0.6 times the expected ratio. This /i/e 
anomaly was interpreted Q as evidence for neutrino oscillations with large amplitude and 
neutrino mass-squared difference 5m 2 atm ~ 10~ 2 eV 2 . Continued experimental studies [[j],0, 
especially by the SuperKamiokande (SuperK) collaboration, have firmly established that the 
deviation of the /i/e ratio from expectation is due to a deficit of muon events. This muon 
deficit increases with zenith angle, and hence with path length, and is consistent with expec- 
tations for muon-neutrino oscillations to some other neutrino flavor or flavors [y^ — > u x ) with 
maximal or near-maximal amplitude and 5m 2 &tra ~ 3.5 x 10~ 3 eV 2 . In principle u x could be 
v e (electron-neutrino), v T (tau- neutrino), or v s (sterile neutrino) ||. However, the observed 
v e flux is in approximate agreement with the predicted v e flux for all zenith angles |l| , which 
rules out z/ M — > v e oscillations with large amplitude. The null results from the CHOOZ and 
Palo Verde reactor v e disappearance experiments also exclude v e — > oscillations at a 
mass-squared-difference scale > eV 2 with amplitude > 0.1. Furthermore, large ampli- 
tude — > v s oscillations at the 5m 2 tm scale are also excluded by SuperK. This is because 
^ii — > v s oscillations are expected to be significantly affected by propagation through mat- 
ter 0H, causing a distortion in the zenith-angle distribution at large angles (corresponding 
to long path lengths) that is not present in the data ||. The zenith-angle distribution ob- 
served by SuperK excludes — > v s oscillations of maximal amplitude at 99% confidence 
level ||. We conclude that, if the oscillation interpretation of the atmospheric neutrino 
deficit is correct, the dominant mode must be — > v T oscillation, with the possibility of 
some smaller amplitude muon-neutrino oscillations to sterile and/or electron- neutrinos ||T0"| . 



An exotic alternative interpretation |]TT| of the atmospheric neutrino disappearance re- 
sults is that a neutrino mass-eigenstate, which is a dominant component of the state, 



decays to a lighter mass-eigenstate and a Majoron [12]. The first oscillation minimum in 



^ ~~ * v v must be observed or excluded to differentiate neutrino oscillations from neutrino 
decays. Unfortunately, the SuperK neutrino-energy and angular-resolution functions smear 
out the characteristic event rate dip that would correspond to the first oscillation mini- 
mum, which cannot therefore be resolved. 

Progress in establishing neutrino oscillations at the atmospheric scale is expected in the 
near future at accelerator neutrino sources with detectors at medium to long baselines. The 
K2K experiment, from KEK to SuperK ||13| , with a baseline L = 250 km and average 
neutrino energy (E u ) = 1.4 GeV, is underway. Their preliminary results are in excellent 
agreement with the oscillation expectations (27 events are observed, whereas 27 events 
would be expected with oscillations and 40 events for no oscillations |]T3[). The MINOS 
experiment from Fermilab to Soudan |l4j] with L = 730 km and (E u ) = 3.5 GeV, which 
begins operation in 2003, is expected to resolve the first oscillation minimum in <-» 
and to search for <-> v e oscillations at the <5m 2 tm scale with an amplitude sensitivity of 



10 . Beginning in 2005, similar physics measurements will be made by the ICARUS fl5 



and OPERA [[RJ experiments with neutrinos of average energy (E u ) ~ 20 GeV from CERN 
detected at L ~ 730 km in the Gran Sasso laboratory. 
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In addition to the atmospheric neutrino deficit, there are other possible indications of neu- 
trino oscillations. In particular, the long-standing deficit of solar neutrinos |17-p0| compared 



to the Standard Solar Model (SSM) predictions |2T] is widely interpreted as an oscillation 
depletion of the v e flux. Note that helioseismology and other solar observations stringently 
limit uncertainties in the central temperature of the sun and other solar model parameters 



22]. The v e deficit relative to prediction is about one-half for the water Cherenkov [18.1 



and Gallium experiments [EO], with the Chlorine experiment WA finding a suppression of 



about one-third. The latest solar neutrino results from SuperK show an electron recoil spec 
trum that is flat in energy, and exhibits no significant day-night or seasonal variation [E3 



An industry has developed to extract the allowed ranges of Sm 2 olllT and v e mixing angles 
that can account for the solar neutrino data. The analyses take account of the coherent 
scattering of v e on matter 0, both in the Sun (the MSW effect [pl| ) and in the Earth ||25| . 



These matter effects can make significant modifications to vacuum oscillation amplitudes. 
Until recently, four viable regions of the parameter space were found in global fits to the 
data: 



(i) LMA — large mixing angle with small matter effects (#m 2 olar rs 10 -5 to 10 _4 eV 

(ii) SMA — small mixing angle with large matter effects {Sm 2 ol&r . ~ 10~ 5 eV 2 ); 

(iii) LOW — large-angle mixing with quasi- vacuum amplitude (Sm 2 olar ~ 10 _7 eV 2 ); 

(iv) VO — large-angle vacuum mixing with small matter effects {5m 2 oXax ~ 10 -10 eV 2 



The latest global solar neutrino analysis by the SuperK collaboration [23[ strongly favors 
solar v e oscillations to active neutrinos [y^ and/or v T ) in the LMA region. A very small area 
in the LOW region is also allowed at 99% C.L., while the SMA and VO regions are rejected 
at 95% C.L. However, other global analyses disagree that the latter regions are excluded ||26|| . 
Moreover, v e — > v s solar oscillations may also still be viable [p7] , p8|| . The relative weighting 
of different experiments (e.g. inclusion or exclusion of the CI data) and whether the 8 B flux 
is held fixed at its SSM value or allowed to float in the global fits presumably account in part 
for the differences in conclusions. It is expected that the KamLAND reactor experiment |29| 



will be able to measure l^mU to ±10% accuracy and sin 2 26*12 to ±0.1 accuracy |JD| if the 
LMA solar solution is correct. 

Finally, the LSND accelerator experiment ]3T[ reports evidence for possible —>■ v e and 
v e oscillations with very small amplitude and 5^ 2 snd ~ 1 eV 2 - If, in addition to the 
LSND observations, the atmospheric and solar effects are also to be explained by oscillations, 
then three distinct 5m 2 scales are needed, requiring a sterile neutrino in addition to the three 
active neutrino flavors [^,|T^,^,^]. The MiniBooNE experiment at Fermilab [32] is designed 
to cover the full region of oscillation parameters indicated by LSND. 

The principal goals of our analyses are to examine the relative merits of different neutrino 
"superbeam" scenarios, where we define a neutrino superbeam as a conventional neutrino 
beam generated by 7r ± decays, but using a very intense megawatt (MW)-scale proton source. 
In particular, we are interested in the physics reach in medium- and long-baseline exper- 
iments with a neutrino superbeam, and how this reach depends upon the beam energy, 
baseline, and the parameters of the neutrino detector. As representative examples we ex- 
plicitly consider water Cherenkov, liquid argon and iron scintillator detectors. However, 
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we note that there is room for new detector ideas, detector optimization, and possibly an 
associated detector R&D program. Therefore, results are presented that apply to any de- 
tector for which the effective fiducial mass and background rates can be specified. Finally, 
we will discuss the role that neutrino superbeams might play |3^j3^] en route to a neutrino 
factory |35],|36|. Our calculations are performed within a three- neutrino oscillation frame- 
work with the parameters chosen to describe the atmospheric neutrino deficit and the solar 
neutrino deficit assuming the LMA solution. However, our considerations for long-baseline 
experiments are relevant even if there are additional short-baseline oscillation effects asso- 
ciated with a fourth neutrino |3"7fl . 

The central objective of long-baseline neutrino oscillation experiments is to determine 
the parameters of the neutrino mixing matrix and the magnitudes and signs of the neutrino 



mass-squared differences; the signs fix the hierarchy of the neutrino mass eigenstates p8 



ye, 



v T ) the mixing matrix relevant to oscillation phenomena can be 



For three neutrinos 

specified by three angles (6*23, 6*12, 6> 13 ) and a phase S associated with CP- violation [see Eq. ([|) 
below]. There are only two independent mass-squared differences (e.g. 5m\ 2 and 8m\^) for 
three neutrinos. The muon-disappearance measurements at SuperK constrain #23 ~ vr/2 
3 x 10~ 3 eV 2 . The other parameter that enters at the leading 8m\ 2 oscillation 



and 5m\ 2 



scale is #13, and its measurement requires the observation of neutrino appearance in v e 



A'" 



Vp — * v K , or v P , — > v T oscillations. 



In long-baseline experiments, if #13 is nonzero, matter effects |f7l, |3"9l , |4"0l modify the proba- 
bility for oscillations involving a v e or v e in a way that can be used to determine the sign of 
f3^.|4lT-44"3"l- Matter effects give apparent CP- violation, but this may be disentangled 

21 at optimally chosen baselines. Matter can also 
Intrinsic CP violation may also be 



Sm 2 32 

from intrinsic CP- violation effects 138 



modify the effects of intrinsic CP or T violation [53 



studied at short baselines where matter effects are relatively small [0,0. CP- violating 
effects enter only for values of L/E u where oscillations associated with the sub leading Sm^i 
become significant [5BLI57]. The most challenging goal of accelerator-based neutrino oscilla- 
tion experiments is to detect, or place stringent limits on, CP violation in the lepton sector. 
We will address the extent to which this may be possible with conventional superbeams. 



II. THREE-NEUTRINO FORMALISM 

The flavor eigenstates v a [a = e, //, r) are related to the mass eigenstates Uj (j = 1, 2, 3) 
in vacuum by 



v, 



■ a = Y. u ^^ (!) 



where U is a unitary 3x3 mixing matrix. The propagation of neutrinos through matter is 
described by the evolution equation |f7lj5"8 



where x = ct and A/2E U is the amplitude for coherent forward charged-current v e scattering 
on electrons, 
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A = 2V2G F Y e pE u = 1.52 x 10" 4 eV 2 x F e p(g/cm 3 ) x E v (GeV) , 



(3) 



where Y e (x) is the electron fraction and p(x) is the matter density. In the Earth's crust the 
average density is typically 3-4 gm/cm 3 and Y e ~ 0.5. The propagation equations can be 
re-expressed in terms of mass-squared differences: 

^ = E ( Sm li u ^ u h + Sm 2 21 U a2 u; 2 + A5 ae 5 Pe ) , (4) 

where Sm 2 k = m 2 — mf,. We assume |5m|]_| <C |5m| 2 |, and that the sign of 5m\ 2 can be 
either positive or negative, corresponding to the case where the most widely separated mass 
eignstate is either above or below, respectively, the other two mass eigenstates. Thus the 
sign of Sm 2 2 determines the ordering of the neutrino masses. The evolution equations can 
be solved numerically taking into account the dependence of the density on depth using the 
density profile from the preliminary reference earth model We integrate the equations 
numerically along the neutrino path using a Runge-Kutta method. The step size at each 
point along the path is taken to be 1% of the shortest oscillation wavelength given by the 
two scales 8m\ 2 and A. 

It is instructive to examine analytic expressions for the vacuum probabilities. We intro- 
duce the notation 

A jk = 5m 2 k L/AE u = 1.27{5m 2 jk /eV 2 {L/km){GeV/E u ) . (5) 

The vacuum probabilities are then given by 

P(y a v p ) = -4:-Re{U a2 U* 3 U; 2 Up 3 ) sin 2 A 32 - ^{U^U^U^U^) sin 2 A 3i 

-^(UaxU^lfaUp) sin 2 A 21 ± 2 JS , (6) 



where J is the CP- violating invariant [60,61 



J = lm(U e2 U: 3 U; 2 U, 3 ) , (7) 

and S is the associated dependence on L and E u , 

S = sin 2 A 2 i + sin 2A 32 - sin 2 A 31 . (8) 

The mixing matrix can be specified by 3 mixing angles (^32,^12)^13) an d a CP-violating 
phase (6). We adopt the parameterization 

/ C13C12 ci 3 si 2 s 13 e~' tS \ 

U = — c 23 Si 2 — s 13 s 23 ci 2 e i5 c 23 Ci 2 — s 13 s 23 Si 2 e* 5 Ci 3 s 23 , (9) 

I iS iS I 

\ s 23 S\ 2 — Si 3 c 23 Ci 2 e — s 23 Ci 2 — Si 3 c 23 Si 2 e Ci 3 c 23 j 

where Cjk = cos9jk and Sjk = sin Ojk- We can restrict the angles to the first quadrant, 
< dij < tt/4, with 5 in the range —71 < 5 < ir. In this parameterization J is given by 

J = si 3 c 2 3 si 2 ci 2 s 23 c 23 s 5 = - sin 29 23 sin 29 12 sin 26 l i 3 ci 3 s (S . (10) 
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For convenience we also define 



K = s 13 c\ 3 s l2 c l2 s 23 c 23 c 6 = - sin 26* 23 sin 26»i2 sin 26 13 c 13 c s . 
Then the vacuum appearance probabilities are given by 



(11) 



PU 



sLs\ 2 sin 2 26»i3 - sin 2 A 32 + s\ndL sin 2 26>i 3 + AK 



sin A 



31 



+ 



c 2 3 (c 2 3 - s 2 3 S2 3 ) sin 2 26>i2 + AK cos 26 12 sin 2 A 21 + 2JS 



C 13l C 12 



s l2 s 2 3 ) sin 2 26*23 + AK cos 6* 



7 23 



sin A 



32 



+ 

+ 2 



c i3( s i2 _ c i2 s i3) sin2 2 ^23 - AK cos 6* 23 sin 2 A 3i 



sin 2 26*0 



23 1*13 



s{ 2 c{ 2 (l + s{ 3 Y + sf 3 sin 2 26 12 (l + sin 2 2 9 23 s 



+ sin 26*12 cos 6*i 2 sin 6* 23 cos 6* 23 Si 3 C5(l + s 2 3 ) sin 2 A 2 i — 2 JS 



c 2 23 s\ 2 sin 2 26*13 + 



sin A 32 + 



C2 3 c 2 2 sin 2 26*13 - 



+ 



'13 



s 2 3 c 2 3 ) sin 2 26*12 - AK cos 26*i 2 



sin 2 A 



21 



sin 2 A 3 i 
2JS 



(12) 



(13) 
(14) 



The corresponding probabilities P{y a — > vp) can be obtained by reversing the sign of 5 in 
the above formulas (only the JS term changes sign in each case). The probabilities for 
v a are the same as those for u a — > v@, assuming CPT invariance. Tests of CPT 
non-invariance are important |62|-|S^| but beyond the scope of the present analysis. The Ajj 
are not independent, and can be expressed in terms of A atm = A 32 and A sun = A 2 i. Then 
A 3i = A atm + A sun and 



sin 2 A 3 i = sin 2 A atm + sin 2 A snn cos 2A aim + - sin 2A sun sin 2A atm , (15) 
S = 2 (sin 2A sun sin 2 A atm + sin 2A aim sin 2 A snn ) . (16) 

Since Eqs. (p~2]) — (|16D in their exact form are somewhat impenetrable, we make a few 
simplifying assumptions to illustrate their typical consequences. First, it is advantageous 
in long-baseline experiments to operate at an L/E u value such that the leading oscillation 
is nearly maximal, i.e. A atm ~ ir/2. Since 8m 2 sun £^ 2 tm , A sun 1 and to a good 
approximation we can ignore terms involving sin 2 A sun . Also, since 6*i 3 is already constrained 
by experiment to be small, for the terms involving A sun we retain only the leading terms 
in 9\ 3 . Second, at the value A atm ~ n/2 for which the leading oscillation is best measured, 
sin2A atm ~ 0. Even if A atm is not close to 7r/2 for all neutrino energies in the beam, an 
averaging over the energy spectrum will suppress sin2A aim if the neutrinos at the middle 
of the spectrum have A atm ~ tt/2. With the above approximations the vacuum oscillation 
probabilities simplify to 



P("» 
P{v„ 



PU 



' e ) ~ sin 2 A atm (s\ 3 sin 2 26* i3 + 4 J sin 2 A 
v) ~ sin 2 A atm (sin 2 26* 23 - 4Jsin2A s , 
v) ~ sin 2 A atm (c 2 3 sin 2 26 13 + AJ sin 2A sm ) 



(17) 
(18) 
(19) 
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It is interesting to compare the relative sizes of the leading CP- violating (CPV) and CP- 
conserving (CPC) terms in the — > v e oscillation probability: 



CPV 4Jsin2A sun / sin 28\2 sin 26*23 \ /A sun sin<5" 



CPC s 2 23 sm 2 28 13 V 2s 2 23 I V 8 



18 



(20) 



For the standard three-neutrino solution to the solar and atmospheric data with large-angle 
mixing in the solar sector, the first fraction on the right-hand side of Eq. (|20|) is of order 
unity, and the relative size of the CPV term is 

CPV A sun sin 5 

CPC 8 13 [ ' 

As an example, with 5m 2 un ~lx 10~ 4 eV 2 and L/E v ~ 300 km/GeV, A sun ~ 0.04; then 
with 5 = vr/2 and sin 2 2# 13 = 0.1 (its maximum allowed value), the CPV term is about 25% 
of the CPC term. Smaller values of Sm^i or sin 5 decrease the ratio in Eq. (pl|); smaller 
values of #13 increase it. 

While smaller values of #i 3 give a larger relative CPV term, they will also reduce the 
overall — > v e event rate since the CPC term is proportional to sin 2 26*13 . The CPV effect 
may be hard to measure if the event rate is low (due to insufficient flux or a small detector). 
Because the number of CPC events is proportional to sin 2 2# 13 , the statistical uncertainty 
on the CPC event rate is proportional to #13 for small #13 and Gaussian statistics. Since 
the number of CPV events is also proportional to #13, the size of the CPV signal relative 
to the statistical uncertainties does not decrease as 8\ 3 becomes smaller. Therefore, a priori 
it does not follow that small #13 automatically makes CPV undetectable p5 . 



On the other hand, even if the event rate is high enough to overcome the statistical 
uncertainties on the signal, backgrounds will limit the ability to measure CPV. Background 
considerations place an effective lower bound on the values of sin 2 26?i3 for which a CPV 
search can be made. This can be quantified by noting that the ratio of the number of CPV 
events N C pv to the uncertainty due to the background is N C pvl V Ib^o (assuming Gaussian 
statistics), where N is the number of events without oscillations and fp is the background 
fraction. Nqpv can be expressed as the CPV part of the oscillation probability times N Q . 
Using the expression for the probability in Eq. (|i~7j), it follows that a 3a CPV effect in the 
Vfj, — > v e channel requires 

sm 2 20 13 > . 2 9 . 2A ^ • (22) 
A Ln sin 5 N 

For 5 = 7r/2 and A sun = 0.03, a typical experiment with fp = 0.01 and iVo = 10 4 can detect 
CPV for sin 2 26*13 > 0.01. The detailed calculations in Sec. confirm this approximate 
result. 

The preceding discussion applies only when the corrections due to matter are not large, 



generally when L is small compared to the Earth's radius. Reference ]54]] gives approximate 
expressions for the probabilities when the matter corrections are small but not negligible. 
However, the most striking matter effects occur when the matter corrections are large and 



the expansions of Ref. [S3| are no longer valid (see, e.g., the plots of oscillation probabilities 



in matter given in Ref. [66] ) 
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Some of the qualitative properties of neutrino oscillations in matter can be determined 
by considering only the leading oscillation and assuming a constant density. There is an 
effective mixing angle in matter defined by 



sin 2 29™ 



sin 2 26 



13 



(^-cos2fl 13 ) + sin 2 29 13 



(23) 



where A is given in Eq. @. The oscillation probabilities in the leading oscillation approxi- 
mation for constant density are pl| , |67 



P{v P 



s\ 3 sin 2 26™ sin 2 A™ , 



sin 6™) 2 sin" A m + (cos 6™Y sin z A™ - (sin 0™ cos 6™Y sin" A 



'13; 



sin 2 26 23 



c 2 3 sin 2 20™ sin 2 A™, 



TTTl \ 2 



m\2 . 



l 21 



'13 



31 



'13 1 



'13; 



\32 



(24) 



where the oscillation arguments are 



Am 
^32 



and 



*atm 



s 



\rn _ a 
^31 — '-^■atm 



A 



+ s 



Am 
^21 



A 



atm ', 



1 + 



.4 



3 m atm 



-S 



(25) 



S 



\ 



A 



cos 26 



13 



+ sin 2 26 13 



(26) 



The A^ term in Piy^ 



u T ) must be retained here because it is not necessarily negligible 
compared to A^, due to matter effects. The expressions for antineutrinos may be generated 
by changing the sign of A. 



In Eq. 



5m atrn COS 20 



13 



3|) there 
(A ~ - 



is a 



resonant enhancement of — > v e oscillations when A ~ 
-Sm 2 tm cos 2# 13 for antineutrinos). This occurs for neutrinos when 



5m 2 atm > and for antineutrinos when 5m 2 atm < 0. On resonance, there is a suppression for 
antineutrinos (neutrinos) when 5m 2 atm > (5m 2 atm < 0). This enhancement of one channel 
and suppression of the other then gives a fake CP violation due to matter effects. 

In the event that the contribution of the sub-leading oscillation is not negligible, the true 
CPV effects due to 5 also enter, but they may be masked by matter effects. Numerical 
calculations [52] show that for distances larger than 2000 km matter effects dominate the 
true CPV for sin 2 26 13 > 0.001 and vice versa for sin 2 26 13 < 0.001. 

As long as sin 2 26\ 3 is not too small, one approach is to have L large enough so that the 
dominant CPV effect is from matter and the sign of Sm 2 tm is clearly determinable; then 
the true CPV effect can be extracted by considering deviations from the CP-conserving 
predictions fj2| . With large L, the neutrino energies must be high enough that A atm ~ n/2 
(e.g., L ~ 3000 km requires E u ~ 10 GeV). An alternative approach is to have short L where 
the matter effects are relatively small ||54)| ; this usually requires a smaller E v to have A atm 
of order 7r/2 (e.g. 
in this paper. 



L ~ 300 km and E v ~ 1 GeV). We will study both of these possibilities 



S 



For sin 2 26^ < 0.001, the matter effect is similar in size or smaller than the true CPV 
effect, and it may not be possible to distinguish between large intrinsic CPV with very 
small #13 from no intrinsic CPV with a moderate-sized #13. Even in experiments at short 
distances (where the matter effect is small) the number of appearance events may be too 
small relative to the background to have a statistically significant difference between the 
neutrino and antineutrino oscillation probabilities. The existence of intrinsic CP violation 
may be very difficult to determine in this case. 

III. CONVENTIONAL NEUTRINO BEAMS 

Conventional neutrino beams are produced using a pion decay channel. If the pions are 
charge-sign selected so that only positive (negative) particles are within the channel, the 
pion decays n + — > /U + z/^ (n~ — > yTv^) will produce a beam of muon neutrinos (antineu- 
trinos). The beams will also contain small components of v e and V e from kaon and muon 
decays. For a positive beam, the dominant decays that contribute to the v e component are 
K + — > 7T°e + z/ e and fi + — > e + u e u fl . If the pion beam has not been charge-sign selected there 
will also be a contribution from K\ — > 7r ± e T z/ e decays. The u e + V e "contamination" can 
be minimized using beam optics that disfavor decays occurring close to the target [note: 
r(ii' ± ) ~ 0.5r(7r ± )], and choosing a short decay channel to reduce the contribution from 
muon decays. These strategies enhance the flavor purity of the beam, but reduce the beam 
flux. Depending on the beamline design, the resulting v e + V e contamination is typically 
at the few parts in 100 to a few parts in 1000 level. The intrinsic v e component in the 
beam produces a background that must be subtracted in a — > u e oscillation search. Ulti- 
mately, the systematic uncertainty associated with the background subtraction will degrade 
the sensitivity of the oscillation measurement. 

To maximize the neutrino flux in the forward direction it is desirable that the pion beam 
divergence is small within the decay channel. The required radial focusing can be provided 
by a quadrupole channel and/or magnetic horns. The beamline optics (dipoles, horns, and 
quadrupoles) determine the peak pion energy and energy spread within the decay channel, 
and hence determine the neutrino spectrum. If the optics are designed to accept a large pion 
momentum spread the resulting wide band beam (WBB) will contain a large neutrino flux 
with a broad energy spectrum. If the optics are designed to accept a smaller pion momentum 
spread, the resulting narrow band beam (NBB) will have a narrower energy spread, but a 
smaller flux. 

A. Detectors and backgrounds 

We are primarily interested in searching for, and measuring, — > v e and — > v e os- 
cillations. The experimental signature for these oscillation modes is the appearance of an 
energetic electron or positron in a charged-current (CC) event. The electron must be sepa- 
rated from the hadronic remnants produced by the fragmenting nucleon. Backgrounds can 
arise from (i) energetic neutral pions that are produced in neutral-current (NC) interactions, 
and subsequently fake a prompt-electron signature, (ii) energetic neutral pions that are pro- 
duced in CC interactions in which the muon is undetected, and the n° fakes an electron, (iii) 
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charm production and semileptonic decay, (iv) — > i/ r oscillations followed by decay of the 
tau-lepton to an electron. Backgrounds (iii) and (iv) can be suppressed using a low-energy 
neutrino beam. 

Background (i) is potentially the most dangerous since leading 7r° production in NC 
events is not uncommon. Indeed, in a recent study [pS| , |6"5|1 using a low energy beam it has 
been shown that in a water Cherenkov detector (e.g. SuperK) it is difficult to reduce this 
background to a level below (9(3%) of the CC rate. A liquid argon detector is believed to 
provide much better 7r°-electron discrimination, and will perhaps enable the ir° background 
to be reduced to 0(0.1%) of the CC rate ]7TJ. Based on these considerations there are 
two different detector strategies. We can choose a water Cherenkov detector, enabling us 
to maximize the detector mass, and hence the event statistics, but obliging us to tolerate 
a significant background from 7r° production in NC events. Alternatively, we can choose a 
detector technology that highly suppresses the n° background, but this will oblige us to use 
a smaller fiducial mass, and hence lower event statistics. 

For a given choice of beamline design, baseline, and detector parameters, the experi- 
mental Ufj, — > v e and — > v e sensitivities can be calculated. It is useful to define some 
representative scenarios, characterized by (a) the parameters of the primary proton beam 
incident on the pion production target, (b) the data sample size D (kt-years), defined as 
the product of the detector fiducial mass, the efficiency of the signal selection requirements, 
and the number of years of data taking, (c) the background fraction defined as the back- 
ground rate divided by the CC rate for events that survive the signal selection requirements, 
and (d) the fractional systematic uncertainty Cf B / /b on the predicted 

We will consider neutrino superbeams that can be produced with MW-scale proton beams 
at low energy (E u ~ 1 GeV) at the proposed Japan Hadron Facility ( JHF) and at high energy 
(E u > 3 GeV) at laboratories with high-energy proton drivers that might be upgraded to 
produce these superbeams (BNL, CERN, DESY, and Fermilab). For the sake of definiteness, 
in the following we will restrict our considerations to two explicit MW-scale primary proton 
beams. First, we will consider the 0.77 MW beam at the 50 GeV proton synchrotron of the 
proposed JHF, and a 4 MW upgrade which we refer to as SJHF (Superbeam JHF). Second, 
we consider a 1.6 MW proton driver upgrade that is under study at Fermilab. With this 
new proton driver, and modest upgrades to the 120 GeV Fermilab Main Injector (MI), it is 
possible to increase the beam current within the MI by a factor of four, and hence increase 
the intensity of the NuMI beam by a factor of four, which we refer to as SNuMI (Superbeam 
NuMI). Higher beam intensities are precluded by space-charge limitations in the MI | |71| . 
For both the SJHF and SNuMI cases we will assume a run plan in which there is 3 years of 
data taking with a neutrino beam followed by 6 years of data taking with an antineutrino 
beam. In principle the high-energy superbeams could be produced at any of the present 
laboratories with high-energy proton drivers. 

We define three aggressive detector scenarios, which are summarized in Table [J 

Scenario A, which might be realized with a liquid argon detector. We choose a 30 kt 
fiducial mass, which has been considered previously for a neutrino factory detector. 
We assume tight selection requirements are used to suppress the n° background, and 
take the signal efficiency to be 0.5. This will result in D = 45 kt-years for neutrino 
running and 90 kt-years for antineutrino running. We assume that backgrounds from 
7T° events contribute 0.001 to [ff^| , and the u e contamination in the beam contributes 
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0.003 to fs in neutrino running, and 0.005 to fs in antineutrino running. We neglect 
all other backgrounds. Hence, fs = 0.004 (0.006) for neutrino (antineutrino) running. 
Finally, we will assume that we know the background rate with a precision of 10% 

(cTfJfB = 0.1). 

Scenario F, which might be realized with a fine-grain iron sampling calorimeter. We 
choose a 10 kt fiducial mass, which is a factor of 10 larger than the THESEUS de- 



tector | 75fl . Since this fiducial mass is smaller than for the alternative scenarios we 
are considering, we will assume that the selection cuts are not tight, and that the 
selection efficiency is 0.9. This will result in D = 27 kt-years for neutrino running 
and 54 kt-years for antineutrino running. We assume that backgrounds from tt° events 
contribute 0.01 to fs, and the v e contamination in the beam contributes 0.003 to fs in 
neutrino running, and 0.005 to Jb in antineutrino running. We neglect all other back- 
grounds. Hence, fs = 0.013 (0.015) for neutrino (antineutrino) running. Finally, we 
will assume that we know the background rate with a precision of 10% (ct/ s //b — 0.1). 

Scenario W, which might be realized for a low-energy beam with a water Cherenkov detec- 
tor. We choose a 220 kt fiducial mass, a factor of 10 larger than the SuperK detector. 



Guided by the study described in Ref. [68] we will assume the selection requirements 



used to suppress the ir° background result in a signal efficiency of 0.68. This will result 
in D = 450 kt-years for neutrino running and 900 kt-years for antineutrino running. 
We assume that backgrounds from n° events dominate, and set fs = 0.02. Note that 
a detailed detector simulation has obtained fs = 0.03 for a water Cherenkov detector 



at a low energy NBB at the JHF ||68|| . With further optimization the choice fs = 0.02 
might therefore be realizable at low energy, but for higher energy (> 1 GeV) neutrino 
beams the rejection against the tt° background is expected to be much worse. Hence, 
a new detector technology might be required for this scenario to make sense at high 
energies. Finally, we will assume that we know the background rate with a precision 
of 10% (cr fB /f B = 0.1). 

Scenarios A, F, and W are very aggressive, and may or may not be realizable in practice. 
In the following we will explore the oscillation sensitivity as a function of D, f B , Of B j '/b, 
baseline, and neutrino beam energy. Scenario F is clearly inferior to scenario A, which has 
larger D and smaller fg. Therefore, in the following we will not discuss the scenario F physics 
potential in detail, but we will indicate the scenario F physics potential on several relevant 
figures. We will use scenarios A and W extensively to illustrate the physics potential of 
upgraded conventional neutrino beams, and facilitate a discussion of the challenges involved 
in probing small values of sin 2 26*13. 



IV. sin 2 2#i3 REACH 

For a given neutrino beam, baseline, and detector, we wish to calculate the resulting 
sin 2 2#i3 reach, which we define as the value of sin 2 2#i3 that would result in a z/ M — > u e 
or — > v e signal that is 3 standard deviations above the background. In our analysis we 
will take into account the Poisson statistical uncertainties on the numbers of signal and 
background events, and the systematic uncertainty on the background subtraction. Our 
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prescription for determining the sin 2 2#i3 reach is given in the appendix. In the following, 
unless otherwise stated, the sin 2 2# 13 reaches are calculated setting the sub-leading oscillation 
parameters sin 2 2# 12 = 0.8 and Sm^i = 10~ 5 eV 2 . The small Sm^i effectively switches off 
contributions from the sub-leading scale. Larger values of Sm^ can yield contributions to 
the Ufj, — > v e signal, but these contributions are not important unless sin 2 2#i3 is significantly 
less than 10 -3 . 



A. Sensitivity at the Japan Hadron Facility 

The Japan Hadron Facility working group has recently investigated |)8| the — > v e 
oscillation sensitivity attainable at the proposed 0.77 MW 50 GeV proton synchrotron in 
Japan, using a 295 km baseline together with the SuperK detector. In their study they 
considered a variety of low energy ((E u ) ~ 1 GeV) WBB and NBB beamline designs. The 



resulting experimental scenario is similar to the one later considered in Ref. p3|| . The 
conclusions from the study were: (i) With a water Cherenkov detector the sensitivity is 
limited by the 7r° background produced in NC events. To minimize the background a 
NBB must be used, since in a WBB the high energy tail will be the dominant source of 
background events, (ii) With an effective — > u e oscillation amplitude of sin 2 29^ e = 0.05 
at an oscillation scale of 8m\ 2 = 0.003 eV 2 , the best set of selection requirements identified 
in the study yielded fs = 0.03, and a signal:background ratio of 1:1 with 12.3 signal events 
per year in the SuperK detector. A 10% systematic uncertainty on the background rate 
was assumed. If no signal is observed after 5 years of running the expected limit would be 
sin 2 29^ e < 0.01 at 90% C.L., which corresponds to a sin 2 2^13 reach of 0.05. (iii) With a 
detector of mass 20 x SuperK, after 5 years running the resulting limit in the absence of 
a signal would be sin 2 2# Me < 0.003 at 90% C.L., which corresponds to a sin 2 2^13 reach of 
0.01. (iv) The energy distribution of the background events is similar to the corresponding 
distribution for the signal. Hence, the z/ M — > v e search is essentially a counting experiment. 
Using these JHF study results, we find that the appropriate values to use in evaluating 
the sin 2 29 13 reach for the JHF to SuperK experiment are fs = 0.03, (Jf B /f B = 0.1, and 
D = 75 kt-years (SuperK with 5 years exposure and a signal efficiency of 68%). With 
these values, our statistical treatment recovers the 90% C.L. results presented in the JHF 
report In our calculations for SJHF we used the neutrino infraction rates presented in 



Ref. [58 



We can now investigate the dependence of the sin 2 2# 13 reach on the detector parameters, 
and hence try to understand whether a massive water Cherenkov detector is likely to be the 
best option. In Fig. |TJ contours of constant sin 2 2#i3 reach are shown as a function of the 
dataset size D and the background rate fs for 3 years of running at the 0.77 MW JHF beam 
(left-hand plots) and at an upgraded 4 MW SJHF beam (right-hand plots). The lower panels 
show how the sin 2 2# 13 reach varies with Oj B j fg. The contours have a characteristic shape. 
At sufficiently large D the sin 2 2# 13 sensitivity is limited by the systematic uncertainties 
associated with the background subtraction, and the reach does not significantly improve 
with increasing dataset size. The contours are therefore vertical in this region of Fig. [I]. At 
sufficiently small D the sensitivity of the — > v e appearance search is limited by signal 
statistics, and further reductions in fg do not improve the sin 2 2# 13 reach. The contours are 
therefore horizontal in this region of Fig. |I|. The positions in the (fs, -D)-plane corresponding 
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to our three detector scenarios (A, F, and W) are indicated on the figure. For the 0.77 MW 
machine the two scenarios (A and W) both yield reaches in the range sin 2 2# 13 ~ 0.015 to 
0.03. However, the water Cherenkov sensitivity is limited by the systematic uncertainty on 
the substantial 7r° background. Hence, the sin 2 2# 13 reach for scenario W does not improve 
substantially when the accelerator beam is upgraded from 0.77 MW to 4 MW (SJHF). On 
the other hand this upgrade would result in a substantial improvement in the reach obtained 
with scenario A, which is not background limited, and therefore has a reach improving almost 
linearly with D. We conclude that, even with SJHF, it will be difficult to observe a —>■ v e 
signal if sin 2 2#i3 is less than about 0.01. This conclusion is consistent with the JHF study 
group analysis, but is in conflict with the expectations of Ref. [[HJ. On the positive side, 
if sin 2 2^13 is larger than 0.01, a 1 GeV neutrino beam at JHF or SJHF would permit the 
observation of a — > v e signal. Detector scenario W does slightly better for a 0.77 MW 
JHF, while scenario A does slightly better for a 4.0 MW SJHF. 

Finally, we consider whether the sin 2 2# 13 reach at a 1 GeV JHF or SJHF neutrino beam 
can be improved with a different choice of baseline. Contours of constant reach in the 
(L, .D)-plane are shown for scenarios A and W in Fig. |2|. A baseline of 295 km does indeed 
yield the optimal reach for the water Cherenkov scenario. For scenario A, a slightly shorter 
baseline (200 km) would yield a slightly improved reach. 



B. Sensitivities for long baseline experiments 

1. Decay channel length restrictions 

Consider next the sensitivity that can be achieved with longer baselines and higher en- 
ergies. We begin by considering how restrictions on the decay channel length reduce the 
neutrino flux for very long baselines. In a conventional neutrino beamline design it is desir- 
able that the pion decay channel is long enough for most of the pions to decay. However, 
for very-long-baseline experiments the decay channel must point downwards at a steep an- 
gle, and the geology under the accelerator site may impose significant constraints on the 
maximum length of the decay channel. In practice, an upgraded long-baseline conventional 
neutrino beam would be sited at an existing particle physics laboratory having a high-energy 
proton accelerator: Brookhaven or Fermilab in the US, CERN or DESY in Europe, or the 
planned JHF laboratory in Japan. The rock characteristics under the JHF site are expected 
to be determined next year by drilling ||70|| . The site with the deepest viable rock layer in 
the US is Fermilab, which sits above approximately 200 m of good rock. The Brookhaven 



and DESY |74| laboratories sit just above the water table — an impediment that would 
have to be overcome before a high-energy long-baseline beam could be proposed. The depth 
of the good rock (Molasse) under CERN varies between about 200 m and 400 m, depending 
on location |72|| . The impact of these restrictions on the maximum decay channel length is 
shown as a function of the baseline length in Fig. |3| for the Fermilab and CERN sites. The 
resulting fraction of pions that decay within the decay channel is summarized in Table [H] for 
several neutrino beam energies. The channel length calculations were performed assuming 
that (i) the proton accelerator is at a depth of 10 m, (ii) the beam is then bent down to point 
in the appropriate baseline-dependent direction using a magnetic channel with an average 
field of 2 Tesla, and (iii) once pointing in the right direction the proton beam enters a 50 m 
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long targeting and focusing section, after which the decay channel begins. The maximum 
decay channel length then depends upon whether the channel extends all the way to the 
bottom of the usable rock layer, or whether this rock layer must also accommodate a near 
detector. Results for both of these cases are presented in Fig. || and Table ||. In the near- 
detector case the maximum decay channel length has been reduced by 100 m to allow for 
the shielding and detector hall. The pion decay fraction estimates have been made assuming 
that all of the decaying pions have the average pion energy in the channel. 

The decay fractions in Table |J show that the site-dependent depth restrictions will 
result in a significant reduction in the neutrino beam intensities for high-energy long-baseline 
beams. For example, at the Fermilab site there is no room for a near detector if the baseline 
is 9300 km (Fermilab to SuperK). With the medium energy beam and a baseline of 7300 km 
(Fermilab to Gran Sasso) only 17% of the pions decay within the channel. Hence, the channel 
length restrictions would exclude, or at least heavily penalize, the extremely-long-baseline 
ideas proposed by Dick et al. [|75[] . Clearly, decay channel length restrictions must be taken 
into account when comparing choices of baseline and beam energy. 



2. sin 2 2#i3 reach for scenarios A and W 

We are now ready to consider the sin 2 2# 13 reach that can be obtained in a long-baseline 
experiment. In our main discussion we consider v e —>■ appearance with b~m\ 2 > 0; the 
5m\ 2 < case is discussed at the end of this section. Our calculations use the WBB spectra 



and interaction rates presented in the MINOS design report [[14J] for the low-energy (LE) horn 
configuration [E v ~ 3 GeV), the medium-energy (ME) horn configuration [E v ~ 7 GeV), 
and the high-energy (HE) horn configuration (E u ~ 15 GeV). After accounting for the 
decay channel length restrictions arising from a maximum depth requirement of 200 m, the 
neutrino fluxes are assumed to scale with the inverse square of the baseline length. 



The calculated reaches are listed in Tables III and IV for detector scenarios A and W, and 



several baselines: L = 730 km (Fermilab — > Soudan or CERN — > Gran Sasso), L = 2900 km 
(Fermilab LBNL/SLAC), L = 7300 km (Fermilab Gran Sasso), and L = 9300 km 
(Fermilab — > SuperK). Note that the shortest baseline (730 km) has a very limited sin 2 26 13 
reach for all the beams, and the lowest energy beam (LE) has a very limited sin 2 2#i3 reach 
for all baselines. The best reach for detector scenario A is sin 2 26*i3 = 0.003, which is obtained 
with a baseline that is not too long (e.g. 2900 km). The best reach for detector scenario 
W is also sin 2 2# 13 = 0.003, and is obtained with long baselines (e.g. 7300 km or 9300 km) 
which benefit from the enhancement of the oscillation amplitude due to matter effects. The 
reaches for the two longest baselines are about the same since the increase of the matter 
enhancement as L increases is compensated by the decrease in the pion decay fraction due 
to the decay channel length restriction. 



To further illustrate the impact of the decay channel length restrictions on the sin 29 
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reach for long-baseline experiments, in Fig. [| contours of constant reach are shown in the 
(J'b, -D)-plane for L = 7300 km with the channel length restrictions (right-hand plots) and 
without the channel length restrictions (left-hand plots). The scenario W point lies in the 
systematics-dominated (vertical contour) region, and is therefore not significantly affected 
by a reduction in D due to the decay channel length restrictions. However, the scenario A 
point lies between the systematics-limited and statistics-limited regions of the plot, and is 
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significantly affected by the reduction in neutrino flux due to the channel length restrictions. 
Indeed, in scenario A the sin 2 2# 13 reach at the HE beam is degraded from about 0.0015 to 
about 0.004 by the channel length restriction. 

3. Dependence on detector parameters 

We can now explore the dependence of the sin 2 2# 13 reach on the baseline, beam energy, 
and detector parameters. In Fig. [| contours of constant sin 2 2# 13 reach are shown in the 
(/b, -D)-plane for L = 730 km and 2900 km. As already noted, for our detector scenarios 
A and W the best reach is ~ 0.003, obtained with scenario A at L = 2900 km using either 
the ME or HE beams (Figs. |^e and ||f), or with scenario W at L — 7300 km using the same 
beams (Figs. |]e and|4]f). It is interesting to consider what improvements to scenarios A and 
W would be required to obtain a reach of 0.001, for example. This goal can be attained 
by decreasing the background fraction for scenario W to fs ~ 0.004 (or alternatively 
increasing the dataset size D for scenario A by a factor of 10) and using the HE beam at a 
very long baseline. The goal could also be attained by decreasing fs for scenario A by an 
order of magnitude and using the high energy beam and a baseline of 2900 km, for example. 
None of these revised detector scenarios seems practical. An alternative strategy is to try 
to find a detector scenario with a smaller systematic uncertainty on Fig. |^ shows, for 
the ME and HE beams, contours of constant sin 2 2#i 3 reach in the (fs, -D)-plane for several 
different af B /fs, and for baselines of 2900 km, 4000 km, and 7300 km. The scenario W 
sensitivity would benefit if the systematic uncertainty on the background could be reduced, 
but even a factor of five improvement in <Jf B /fB would not permit a reach of 0.001 to be 
attained. 

Since detector scenarios A and W are ambitious, we can ask what happens if D, J'b, 
or o~ j B / f b must be relaxed. The best reaches obtained with scenario W were for the ME 
and HE beams at very long baselines (e.g. L = 7300 km). In these cases the reach is not 
very sensitive to D, but degrades roughly linearly with increasing f B (Fig. [|) or <Jf B /f B 
(Fig. |6|). Hence, if the achievable background rate is really fs = 0.1 then the sin 2 2# 13 reach 
is well above 0.01 for the observation of a — > v e signal at 3 standard deviations above the 
background. The best reaches obtained with scenario A were for the ME and HE beams at 
long baselines (e.g. L = 2900 km). In these cases the reach is sensitive to both decreases in 
D and increases in The reach can be degraded by a factor of 2 by either reducing D by 
about a factor of 4, or by increasing f B by about a factor of 3 (Fig. g). 

So far we have considered only a few discrete baseline lengths. To explore the reach 
that can be obtained with other baseline choices, Fig. [j] shows, for each of the three NuMI 
beam energies, contours of constant sin 2 2#i 3 reach in the (L, D)-plane for scenarios A and 
W. For scenario A, where backgrounds are less important, the optimal distance varies with 
beam energy; crudely speaking, the optimal L is given by making the vacuum oscillation 
argument 1.27<5m 2 2 L/ (E u ) of order tt/2. For scenario W the backgrounds are more important 
and larger distances give a better sin 2 2#i 3 reach for all three upgraded NuMI beams. 

Finally, we have also studied neutrino beams with higher energy than NuMI. For example, 
the CNGS beam [I7?J at CERN has an average neutrino energy of about 20 GeV. We find that 
for the expected 3 x 10 19 protons on target per year, three years of running will at best give 
a sin 2 26* 13 reach of about 0.01 for either scenario A or W. Upgrading the proton intensity 
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by a factor of four improves the sin 2 2^13 reach to about 0.005. Therefore we conclude that 
the higher-energy CNGS superbeams have similar capability to the SNuMI beams. 

4- Summary o/sin 2 29i% reaches for Sm^ 2 > and 8m\ 2 < 

In summary, Figs. |], [|, and [7] show that the best sin 2 2613 reach that can be obtained 
with detector scenarios A and W is about 0.003. This optimum reach can be obtained in 
scenario A with L ~ 2000-4000 km for the NuMI ME beam or with L ~ 3000-6000 km 
for the HE beam, or in scenario W with L ~ 7000-9000 km for either the ME or HE 
beams. Scenarios A and W require ambitious detector parameters. To improve the reach to 
0.001, for example, requires substantial improvements in <jf B /B, and/or D, and does not 
therefore seem practical. If the scenario A and W parameters cannot be realized the reach 
will be degraded. In particular, a significant increase of fs (or o~f B j js) in either scenario A 
or W would result in a significant decrease in sin 2 2#i3 reach. A significant decrease in the 
data-sample size in scenario A will also degrade the sin 2 26*13 reach. 

Up to now we have considered the sensitivity of long baseline experiments if 8m\ 2 > 0. 
We now turn our attention to the alternative case: 8m\ 2 < 0. In this case long baseline 
experiments using a neutrino beam will suffer from a suppression of signal due to matter 
effects. Therefore, in our scenarios A and W, if no signal is observed after 3 years of neutrino 
running the beam is switched to antineutrinos for a further 6 years of data taking. For 
antineutrino running with 5m| 2 < the results shown in Figs. |], [5], and [7] must be modified 
since the antineutrino cross section is about half of the neutrino cross section. Hence we 
must double the required values on the .D-axes in the various figures. Other modifications 
to the contour plots for antineutrino running with 5m\ 2 < are minor since the matter 
enhancement in this case is similar to the enhancement for neutrinos when 8m\ 2 > (they 
are the same in the limit that the sub- leading oscillation can be ignored). However, the 
positions of the scenario A and W points on the various figures must be moved to account 
for the larger values of D and fg (and potentially er/ B / js). Note that the larger background 
rate associated with antineutrino running in Scenario A will degrade the ultimate sin 2 29 13 
reach for 8m\ 2 < 0; the best reach becomes ~ 0.004. 

V. NEUTRINO MASS HIERARCHY AND CP- VIOLATION 

In the 1 GeV and multi-GeV superbeam scenarios that we have considered it will be dif- 
ficult to observe a — > v e or — > v e signal if sin 2 2#i3 is smaller than about 0.01. However, 
if sin 2 2#i3 is (9(0.01) a z/ M — > v e signal would be observable provided a sufficiently massive 
detector with sufficiently small background is practical. We would like to know if, in this 
case, the sign of 8m\ 2 can be determined in the long-baseline multi-GeV beam experiment, 
and whether CP violation might be observed in either the long-baseline multi-GeV beam 
experiment or the 1 GeV intermediate baseline experiment. We begin by considering the 
CP sensitivity at the SJHF, and then consider the sensitivity for determining CP violation 
and/or the pattern of neutrino masses at long baselines. 
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A. CP violation with a JHF superbeam 



In our SJHF scenario, a CP violation search would consist of running for 3 years with 
a neutrino beam and measuring the number of — ► v e signal events [N(e— )], and then 
running for 6 years with an antineutrino beam and measuring the number of — ► v e 
signal events [iV(e+)]. In our calculations we assume that the antineutrino cross section 
is about one-half of the neutrino cross section, and that the antineutrino flux is the same 
as the neutrino flux. In the absence of CP violation (5 = or 180°), after correcting 
for cross-section and flux differences, we would therefore expect N(e+) ~ N(e—). In the 
presence of maximal CP violation with 5 = 90° [—90°] we would expect N(e+) > N(e— ) 
[N(e+) < N(e— )]. The magnitude of the deviation from N(e+) = N(e—) induced by 
CP violation is quite sensitive to the sub-leading scale dm^. Setting 5 = 90°, in Fig. |8| 
the predicted positions in the [N(e— ), A(e+)]-plane are shown for scenarios A (left-hand 
plots) and W (right-hand plots) The predictions are shown as a function of both sin 2 29 13 and 
5m\ v The error ellipses around each point indicate the measurement precision at 3 standard 
deviations, taking into account both statistical and systematic uncertainties, and using the 
statistical prescription described in the appendix. An overall normalization uncertainty 
(which could account for uncertainties in the flux and/or cross sections) of 2% is included, 
although its effects are generally small. CP violation can be established at the 3a level if 
the error ellipses do not overlap the CP-conserving curves (solid lines, 5 = 0). The curves 
for the other CP-conserving case (5 = 180°) lie very close to the 6 = curves and are not 
shown. 

Note that for scenario W with the upgraded 4 MW SJHF beam, if sin 2 29 13 = 0.1 
(larger values are already excluded), Srn^ = 5 x 10~ 5 eV 2 , and 5 = 90°, then the predicted 
point in the [A(e— ), A(e+)]-plane is just 3cr away from the CP conserving (N(e+) = 
N(e—)) prediction. Alternatively, if sin 2 2# 13 = 0.02, 5rn^ 1 = 1 x 10~ 4 eV 2 (larger values 
are improbable), and 5 = 90°, then the predicted point is also just 3a away from the 
CP-conserving prediction. Hence there is a small region of the allowed parameter space 
(sin 2 2#i3 > 0.02 and <5m 21 > 5 x 10 -5 eV 2 ) within which maximal CP violation might be 
observable at an upgraded JHF if 5m| 2 is in the center of the presently favored SuperK region 
and sin 2 2# 23 ~ 1- It is also possible to detect maximal CP violation for sin 2 2# 13 > 0.05 and 
5rri2i > 10~ 4 eV 2 with the 0.77 MW JHF in the W scenario. Generally detector scenario W 
does better for CP violation, except scenario A is slightly better for sin 2 2#i3 ~ 0.02 at the 
4 MW SJHF. Because the matter effect is small at L = 295 km, predictions for 5m\ 2 > 
and 5m\ 2 < are nearly the same, and hence the sign of 5m| 2 cannot be determined. 

B. CP violation and the sign of 5m 2 2 at long-baseline experiments 

Consider next long-baseline experiments using multi-GeV neutrino beams. The approx- 
imate equality N(e+) ~ N(e—) will be modified by intrinsic CP violation and by matter 
effects. Predictions in the [N(e— ), A(e+)]-plane are shown in Fig. ^ for scenario A using the 
1.6 MW LE superbeam for two values of Sm^i (5 x 10~ 5 eV 2 and 1 x 10~ 4 eV 2 ), and for two 
baselines [L = 730 km and 1800 km). The predictions for each of these cases are shown as a 
function of sin 2 2#i 3 , 5, and the sign of <5m| 2 , with |5m| 2 | = 3.5 x 10~ 3 eV 2 and sin 2 26* 2 3 = 1. 
Note that at L = 730 km the magnitude of the modifications of the appearance rates due 
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to matter effects are comparable to the magnitudes of the modifications due to maximal 
intrinsic CP violation. Furthermore, the expected precisions of the measurements, shown 
on the figure by the 3a error ellipses, are also comparable to the sizes of the predicted CP 
and matter effects. 

Matter effects will cause the two CP-conserving cases 8 = and 8 = 180° to give 
different predictions for iV(e + ) and N(e~), and therefore to establish CP violation the 
signal must be distinguishable from both 5 = and 8 = 180°. Hence, in the scenario we are 
considering, Fig. |5] shows that superbeam measurements with the LE beam at 730 km can 
help to constrain the parameter space, but generally cannot provide unambiguous evidence 
for intrinsic CP violation, and cannot unambiguously determine the sign of 8m\ 2 . The 
only exception to this is if 8m\ 2 > and 8 = —90° (or 8m\ 2 < and 8 = 90°), in which 
case CP violation could be established and the sign of 8m\ 2 determined for sin 2 2#i3 > 
0.02. The CP and matter effects are better separated at L = 1800 km, for which an 
unambiguous determination of the sign of 8m\ 2 seems possible provided sin 2 2# 13 > 0.02, 
although CP violation cannot be established for 8m 21 < 10~ 4 eV 2 . At smaller values of 
sin 2 2#i3 modifications to the appearance rates cannot distinguish between matter and CP 
effects. Note that, because of the matter effect, at distances longer than 1000 km the values 
of 8 that give the largest disparity of N(e + ) and N(e~) are no longer ±90°. Also note 
that the sign of 8m\ 2 is most easily determined when the CPV and matter effects interfere 
constructively to give a greater disparity of N(e + ) and N(e~), and more difficult when 
the CPV and matter effects interfere destructively [i.e., N(e + ) and N(e~) are more equal]. 
Going to even longer baselines, predictions in the [N(e— ), N(e+)] plane are shown in Fig. ^ 
for scenario A (left-hand plots) and scenario W (right-hand plots) with L = 2900 km. The 
predictions are shown for the LE beam (top plots), ME beam (middle plots), and HE beam 
(bottom plots). In general, the sign of 8m\ 2 can be determined provided sin 2 29 13 > 0.02, 
but in none of the explored long-baseline scenarios can CP- violation be unambiguously 
established for 8m^ < 10 -4 eV 2 . 



C. CP- violation and the sign of 5m| 2 at a neutrino factory 

We can ask, how do the CPV and 5m| 2 -sign capabilities of superbeams compare with 
those of a neutrino factory? The relevant experimental signature at neutrino factory is the 
appearance of a wrong-sign muon indicating v e — > (or u e — > v^) transitions. This is a 
much cleaner signature than electron appearance with a superbeam. Hence, background 
systematics are under better control at a neutrino factory, and the expected error ellipses in 
the [N(/i+), N(/x— )]-plane are therefore much smaller. 

In our analysis we assume a 20 GeV neutrino factory with 1.8 x 10 21 useful /i + decays 
(which might be achieved in three years running at at high-performance neutrino factory) 
and 3.6 x 10 21 useful pT decays, a 50 kt iron-scintillator detector pi at distances L = 



1800 km, 2900 km, and 4000 km from the source. For comparison, the total neutrino flux 
for three years running at a distance of 1 km from the source is 2 x 10 19 /m 2 for the neutrino 
factory scenario, while it is 7 x 10 16 /m 2 for SJHF and 4 x 10 18 /m 2 for the SNuMI HE 
beam. We choose an iron-scintillator detector for the neutrino factory analysis since it is 
particularly well-suited for the detection of muons and can be made larger than, e.g., a 
liquid argon detector, at a similar or lower cost. We also take fs = 10~ 4 , <Jf B j Jb = 0.1, and 
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a normalization uncertainty of 2%. This background level can be achieved with a 4 GeV cut 
on the detected muon, which gives a detection efficiency of about 73%, implying an effective 
data sample of D = 110 kt-yr for three years running. 

The corresponding neutrino factory predictions in the [N(/x— ), N (//+)] -plane are shown in 
Fig. [II]. The 1800 km baseline is too short, since matter and CP effects are indistinguishable 
in most cases. At 2900 km the predictions allow an unambiguous determination of the sign 
of 5m^ 2 for much of the parameter space, and the possibility of establishing the existence 
of CPV. At 4000 km the statistical uncertainties are larger, and impair the sensitivity to 
observe CPV. However, matter effects are also larger, and an unambiguous determination 
of the sign of 5m\ 2 is possible down to sin 2 2# 13 of a few x 10~ 4 . For very long baselines (e.g. 
L = 7300 km, Fig. |12|) there is negligible sensitivity to CPV or to 8m\ x , matter effects are 
large, and the sin 2 2#i3 reach for determining the sign of 8m\ 2 approaches 10 -4 . 

VI. SUMMARY 

We have explored the oscillation-physics capabilities of 1 GeV and multi-GeV neutrino 
beams produced at MW-scale proton accelerator facilities (neutrino superbeams). Specif- 
ically, the limiting value of sin 2 2# 13 that would permit the first observation of — > u e 
and/or — + v e oscillations at 3 standard deviations is considered, along with the ability of 
these intense conventional neutrino beams to determine the pattern of neutrino masses (sign 
of 5mfj 2 ) and discover CP- violation in the lepton sector. The figures in this paper provide 
a toolkit for accessing the physics capabilities as a function of the detector specifications, 
characterized by the dataset size D (kt-years) and the uncertainty on the background sub- 
traction (given by f B and af B / '/#). Table |V| summarizes the physics capabilities of some 
beam-detector combinations. Also shown in the table are similar results for an entry-level 
and high-performance neutrino factory with = 20 GeV. 

Determining the optimum detector technology and characteristics is beyond the scope 
of this paper, and may require a detector R&D program. However, for some ambitious but 
plausible detector scenarios we find: 

(i) With a sufficiently ambitious detector, if sin 2 2^13 > few xl0~ 3 and <5m 2 2 > 0, then 

— > v e an d Vfi — > v e signals should be observable at a superbeam. The reach is 
slightly worse if 5m\ 2 < 0. The best reach is obtained with a long-baseline multi-GeV 
superbeam; for example, with the SNuMI ME or HE beams and a baseline > 2000 km. 
This would permit the tightening of constraints on the oscillation parameter space. 
It is important to account for decay channel length restrictions when assessing the 
capabilities of very-long-baseline experiments. 

(ii) If CP is maximally violated in the lepton sector, there is a small region of allowed 
parameter space in which an experiment at a JHF or SJHF beam (E v ~ 1 GeV) might 
be able to establish CP- violation at 3 standard deviations. Except for certain small 
regions in parameter space where matter and CPV effects constructively interfere, a 
long-baseline experiment with conventional superbeams would be unable to unambigu- 
ously establish CP violation because matter effects can confuse the interpretation of 
the measurements. 



19 



(iii) With a sufficiently ambitious detector, if sin 2 29i3 > 0(0.01) there is a significant 
region of parameter space over which a long baseline experiment with a multi-GeV 
neutrino superbeam could unambiguously establish the sign of 5mJ 2 . 

(iv) Lower-energy superbeams do best at shorter distances, with a fair reach for — > v e 
appearance and some CPV capability, but little or no sensitivity to the sign of 5m| 2 ; 
higher-energy superbeams do best at longer distances, with good reach for — > v e 
appearance and sign(5m 2 2 ) determination, but little or no sensitivity to CPV . 

(v) A neutrino factory can deliver between one and two orders of magnitude better reach in 

sin 2 29 13 for v e — > appearance, the sign of Sm 32 , and CP violation; for L ~ 3000 km 
there is excellent sensitivity to all three observables. 

Note that in this study we have restricted our considerations to 1 GeV and multi-GeV 
neutrino beams. The potential of sub-GeV beams is currently under consideration f77||78| . 
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APPENDIX 



To implement the Poisson statistical uncertainties in our analysis of the sin 2 26^ reach 
we use an approximate expression for the upper limit on the number of events from 
the observation of N events, 

Af/- N + SVNTl + (S 2 + 2)/3, (27) 

where S is the number of standard deviations corresponding to the limit. This expression 
gives the correct \jj with an accuracy that is better than 10% for iV < 4, and better than 1% 
for larger N [f79| . If the number of predicted background events is B, the expected number 



of signal events corresponding to an observation 3 statistical standard deviations above the 
background is given by 



N s = 3VB + 1 + 11/3 . (28) 

Let the systematic uncertainty on B be given by U. To account for this systematic uncer- 
tainty, we add it in quadrature with the statistical uncertainty. Defining the quantity 



K = V^ 2 + 9^ 2 , (29) 

the sin 2 2# 13 reach can then be estimated by finding the value of sin 2 26^3 that yields N' s 
signal events. 

To determine the sign of 5m\ 2 and/or search for CP violation with conventional z/ M 
and beams, we will need to compare the — > v e and z/ M — > v e appearance rates (for a 
neutrino factory with a detector that measures muons, we compare the v e — > and v e — > 
appearance rates). As in the case of the sin 2 2# 13 reach, we will be considering the 3a allowed 
regions. 

Let N and N be the number of events that satisfy the signal selection criteria and are 
recorded respectively during neutrino and antineutrino running. If N th and N th are theoret- 
ical predictions for N and N, the region of the N th -N th space allowed by the measurements 
is described by 

(N th -N\ 2 (N th -N\ 2 , s 

+ — <1, (30) 



where aN and a N are the experimental uncertainties on iV and N, respectively. In the 
absence of systematic uncertainties, and in the approximation of Gaussian statistics, the 
3a values are ajy = 3a//V and a N = 3VW. However, since iV and N might be small 
Gaussian statistics may be inappropriate. Instead, we define and a N to correspond 
to the appropriate 99.87% confidence level deviations from the central values of N and N, 
respectively, using Poisson statistics. The expressions for aN and a N will depend on whether 
we are considering an upper or lower limit. 

Consider first the case of an upper limit. We can compute aN using Eq. (|27D, with S = 3, 
yielding: 



aj per = 3VN + 1 + 11/3. (31) 
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To compute the value for for a lower bound, we need an expression for the Poisson lower 
limit given the observation of N events. We use the expression from Ref. [ 791 , namely: 



N (l \- ^= + P NJ ) > ( 32 ) 



9N 3VN 

where with S = 3 we have f3 = 0.222 and 7 = —1.88. This approximate expression for the 
Poisson lower limit on N is accurate to a few percent or better for all N. Hence 

lower = N _ N f 1 _J_ 3 222^-1.88^ (33) 

N \ 9N 3VN J V ; 

The corresponding values for can be found by substituting N for N in Eqs. (|31| ) and 



In practice N and N will contain background components B and B. The predicted 
backgrounds will have associated systematic uncertainties U and U. In this case we can still 
use Eq. (|3Up to determine the allowed regions, but to take account of the background and 
systematic uncertainties the a 2 N and a% are replaced with the substitutions: 

a 2 N -> a 2 N + 9U 2 (34) 



a% -> a% + 9U 2 . (35) 

Other systematic uncertainties on the predicted N and N (for example, the uncertainty on 
the neutrino and antineutrino cross-sections) can be handled in a similar way, by replacing 
otN (ctjv) with the quadrature sum of (ajv) and the additional 99.87% C.L. uncertainty 
on N (N). 
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TABLE I. Parameters for scenarios A, F, and W discussed in the text. The scenarios assume 3 
years of neutrino running and 6 years of antineutrino running. The dataset sizes D and background 
fractions Jb are defined for the event samples after the signal selection requirements have been 
applied. 







Scenario A 


Scenario F Scenario W 






V 


V 


V 


V V 


V 


Fiducial mass 


(kt) 


30 


30 


10 


10 220 


220 


D (kt-years) 




45 


90 


27 


54 450 


900 


Backg. frac. / 


B 


0.004 0.006 


0.013 


0.015 0.02 


0.02 


Backg. uncertainty <t/ s //b 


0.1 


0.1 


0.1 


0.1 0.1 


0.1 


TABLE II. 


The fraction of pions decaying 


5 in a 


channel with the maximum length permitted 


by the depth of viable rock (D max ) under the accelerator site, tabulated as a function of baseline 


L for configurations with and without a near detector. 






SNuMI 


-E^peak) 


-^max 




L 


/decay 


Beam 


(GeV) 


(m) 




(km) 


with near 


no near 


LE 


3 


200 




2900 


0.93 


0.95 










7300 


0.36 


0.56 










9300 




0.37 


ME 


7 


200 




2900 


0.67 


0.72 










7300 


0.17 


0.29 










9300 




0.17 


HE 


15 


200 




2900 


0.41 


0.44 










7300 


0.09 


0.16 










9300 




0.09 


LE 


3 


400 




2900 


0.98 


0.99 










7300 


0.82 


0.88 










9300 


0.67 


0.77 


ME 


7 


400 




2900 


0.93 


0.94 










7300 


0.56 


0.64 










9300 


0.38 


0.47 


HE 


15 


400 




2900 


0.69 


0.72 










7300 


0.29 


0.34 










9300 


0.20 


0.26 
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TABLE III. sin 2 20 is reach (corresponding to a signal that is 3 standard deviations above the 
background after 3 years of running with a neutrino superbeam) shown as a function of baseline 
L for Scenario A described in the text. The oscillation probability — > v e ) corresponding to 

sin 2 2#i3 = 0.01 and the expected numbers of signal events S and background events B are also 
listed. The calculations assume Ara| = 3.5 x 10~ 3 eV 2 , Sm^ = 5 x 10~ 5 eV 2 , and 5 = 0. 



HAT TV TT 

SNuMI 


i^(peak) L 


P 


S 


B 


sin 2 2#i3 reach 


Beam 


(GeV) (km) 










LE 


3 730 


0.0024 


210 


340 


0.006 




2900 


0.0045 


26 


24 


0.008 




7300 


0.012 


4.1 


1.3 


0.02 




9300(*) 


0.016 


3.2 


0.8 


0.02 


ME 


7 730 


0.0016 


370 


910 


0.01 




2900 


0.0075 


120 


62 


0.003 




7300 


0.025 


15 


2.4 


0.006 




9300(*) 


0.035 


14 


1.6 


0.006 


HE 


15 730 


0.0006 


290 


2000 


0.02 




2900 


0.0054 


180 


130 


0.003 




7300 


0.024 


25 


4.2 


0.004 




9300(*) 


0.032 


25 


3.1 


0.004 


(*) No near detector 










TABLE IV. sin 2 29 13 reach (corresponding to a si| 


mal that is 


3 standard deviations above the 


background after 3 years of running with a neutrino superbeam} 


i shown as a function of baseline 


L for Scenario W described in the text. 


The oscillation probability P{v^ — > v e 


) corresponding to 


sin 2 26*13 


= 0.01, and the expected numb 


<ers of signal events S and background events B are also 


listed. The calculations assume Am 2 2 = 


3.5 x 10~ 3 eV 2 , 5m 2 2l = 


5 x 10" 5 eV 2 , 


and 6 = 0. 


SNuMI 


^(peak) L 


P 


S 


B 


sin 2 26*13 reach 


Beam 


(GeV) (km) 










LE 


3 730 


0.0024 


2100 


17000 


0.03 




2900 


0.0045 


260 


1200 


0.02 




7300 


0.012 


41 


67 


0.009 




9300(*) 


0.016 


32 


40 


0.008 


ME 


7 730 


0.0016 


3700 


46000 


0.05 




2900 


0.0075 


1200 


3100 


0.008 




7300 


0.024 


150 


120 


0.003 




9300(*) 


0.035 


140 


80 


0.003 


HE 


15 730 


0.0006 


2900 


98000 


0.1 




2900 


0.0054 


1800 


6700 


0.01 




7300 


0.024 


250 


210 


0.003 




9300(*) 


0.032 


250 


160 


0.003 



(*) No near detector 
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TABLE V. Summary of the sin 2 29\s reach (in units of 1CT 3 ) for various combinations of 
neutrino beam, distance, and detector for (i) a 3a — > v e appearance with Sm^i = 10~ 5 eV 2 , (ii) 
a 3a determination of the sign of Sm^ 2 with Sm^i = 5 x 10 -5 eV 2 , and (hi) a 3a discovery of CP 
violation for three values of Sm^i (in eV 2 ). Dashes in the sign of <5m 2 2 column indicate that the 
sign is not always determinable. Dashes in the CPV columns indicate CPV cannot be established 
for sin 2 2^i3 < 0.1, the current experimental upper limit, for any values of the other parameters. 
The CPV entries are calculated assuming the value of 8 that gives the maximal disparity of N(e + ) 
and N(e~); for other values of 5, CP violation may not be measurable. 



sin 2 2#i3 reach (in units of 10 3 ) 









^ Ve 


Unambiguous 


Possible 3a CPV 








app ear ance 


3a si 


gn((5m§ 2 ) 




Smli (in eV 2 


) 


Beam 


L (km) 


Detector 


8m\ x = 1(T 5 


bm\ x -- 


= 5 x 10" 5 


5 x 10" 


~ 5 1 x 10" 4 


2 x 10~ 4 


JHF 


295 


A 


25 










25 






W 


17 








40 


8 


SJHF 


295 


A 


8 








5 


3 






W 


15 






100 


20 


5 


SNuMI LE 


730 


A 


7 






100 


20 


4 






W 


30 










40 


SNuMI ME 


2900 


A 


3 




6 






100 






W 


8 




15 










7300 


A 


6 




6 












W 


3 




3 








SNuMI HE 


2900 


A 


3 




7 




100 


20 






W 


10 




15 










7300 


A 


4 




4 












W 


3 




3 








20 GeV NuF 


2900 


50 kt 


0.5 




2.5 




2 


1.5 


1.8 x 10 20 fj,+ 


7300 




0.5 




0.3 








20 GeV NuF 


2900 


50 kt 


0.1 




1.2 


0.6 


0.4 


0.6 


1.8 x 10 21 fj+ 


7300 




0.07 




0.1 
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FIG. 1. Contours of constant sin 2 2^13 reach that correspond to a v e —* z/„ signal that is 3 
standard deviations above the background. The contours are shown in the (D, /^)-plane, where 
-D is the data-sample size and Jb the background rate divided by the total CC rate. The contours 
are shown for the 0.77 MW (left-hand plots) and 4.0 MW (right-hand plots) JHF scenarios with 
L = 295 km. The top panels show curves for ot B j /b = 0.1, while the bottom panels show curves 
for (Jf B / fs = 0.1, 0.05, and 0.02. The positions corresponding to the three standard detector 
scenarios defined in Table | are indicated. 
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FIG. 2. Contours of constant sin 2 2#i3 reach that correspond to a v e — > signal that is 3 
standard deviations above the background. The contours are shown in the (D, L)-plane, where 
D is the data-sample size and L the baseline. The panels show predictions for the JHF scenario 
described in the text, for detector scenarios A (left-hand plots) and W (right-hand plots), and for 
0.77 MW (top plots) and 4.0 MW (bottom plots) proton drivers. The positions corresponding to 
scenarios A and W (see Table |) at L = 295 km are indicated. 
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FIG. 3. Maximum length of the pion decay channel that fits within a rock layer that is 200 m 
deep (left-hand plot) and 400 m deep (right-hand plot) shown as a function of baseline. The 
calculation is described in the text. The solid (broken) curves shows the results without (with) a 
near detector. For comparison, the horizontal solid line indicates the NuMI decay channel length. 
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FIG. 4. Contours of constant sin 2 26±3 reach that correspond to a v e — » i/„ signal that is 3 
standard deviations above the background, at L = 7300 km. The contours are shown in the 
(D, /s)-plane, where D is the data-sample size and /b the background rate divided by the total 
CC rate. The contours are shown for the LE (top plots), ME (center plots), and HE (lower 
plots) upgraded SNuMI beams, both with (right plots) and without (left plots) the decay length 
constraint. The systematic uncertainty o"/ s //b = 0.1. The positions corresponding to the three 
standard scenarios defined in Table | are indicated. 



33 




FIG. 5. Contours of constant sin 2 2#i3 reach that correspond to a ^ — ► signal that is 
3 standard deviations above the background, at L = 730 km (left plots) and 2900 km (right 
plots). The contours are shown in the (D, /s)-plane, where D is the data-sample size and Jb the 
background rate divided by the total CC rate. The contours are shown for the LE (top plots), ME 
(center plots), and HE (lower plots) upgraded SNuMI beams. The systematic uncertainty on the 
background subtraction is o"/ s //b = 0.1. The positions of the three standard scenarios defined in 
Table | are indicated. 
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FIG. 6. Contours of constant sin 2 2#i3 reach that correspond to a v e — > f M signal that is 3 
standard deviations above the background, at the upgraded SNuMI ME (left) and HE (right) 
beams. The contours are shown in the (D, /^)-plane, where D is the data-sample size and fs the 
background rate divided by the total CC rate. The contours are shown for L = 2900 (top), 4000 
(center), and 7300 km (bottom). Curves are shown for systematic uncertainties on the background 
subtraction o-/ s //s = 0.1, 0.05, and 0.02. The positions of the three standard scenarios defined in 
Table | are indicated. The decay length constraints have been imposed for L = 4000 and 7300 km. 
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FIG. 7. Contours of constant sin 2 2#i3 reach that correspond to a v e — > signal that is 3 
standard deviations above the background. The contours are shown in the (D, L)-plane, where D 
is the data-sample size and L the baseline. The panels show predictions for the upgraded SNuMI 
LE (top), ME (center), and HE (bottom) beams, and for detector scenarios A (left plots) and 
W (right plots). The decay length constraint is included. The systematic uncertainty on the 
background subtraction <t/ s //b = 0.1. The positions of the standard scenarios defined in Table | 
are shown at L = 730, 2900, and 7300 km. 



36 



A W 




50 100 150 200 250 300 350 500 1000 1500 2000 2500 3000 3500 
N(e-) N(e-) 



FIG. 8. 3<r error ellipses in the [N(e + ), A f (e~)]-plane, shown for the 0.77 MW JHF (top plots) 
and 4 MW SJHF(bottom plots) scenarios with at L = 295 km. The contours are shown for detector 
scenarios A (left) and W (right), with sin 2 26 13 = 0.02, 0.05, and 0.1. The solid (dashed) [dotted] 
curves correspond to 6 = 0° (90°) [—90°] with 8m%x varying from 2 x 10~ 5 eV 2 to 2 x 10~ 4 eV 2 . The 
error ellipses are shown for three simulated data points at = 5 x 10 -5 , 10 -4 and 2 x 10 -4 eV 2 . 
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FIG. 9. 3ct error ellipses in the [N(e + ), iV(e~)]-plane, shown for detector scenario A at the 
upgraded LE SNuMI beam with L = 730 (top plots) and 1800 km (bottom plots). The contours 
are shown for = 5 x 10~ 5 (left) and 10~ 4 eV 2 (right). The solid and long-dashed curves 

correspond to the CP conserving cases 5 = 0° and 180°, and the short-dashed and dotted curves 
correspond to two other cases that give the largest deviation from the CP conserving curves; along 
these curves sin 2 2#i3 varies from 0.001 to 0.1, as indicated. 
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FIG. 10. 3<7 error ellipses in the [N(e + ), iV(e~)]-plane, shown for detector scenarios A (left) 
and W (right) at L = 2900 km with the upgraded LE (top), ME (center), and HE (bottom) 
SNuMI beams. The contours are shown for bm^ = 10~ 4 eV 2 . The solid and long-dashed curves 
correspond to the CP conserving cases 5 = 0° and 180°, and the short-dashed and dotted curves 
correspond to two other cases that give the largest deviation from the CP conserving curves; along 
these curves sin 2 2#i3 varies from 0.001 to 0.1, as indicated. 
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FIG. 11. 3<t error ellipses in the [iV(/i + ), iV(/x )]-plane, shown for a neutrino factory delivering 
3.6 x 10 21 useful decays of 20 GeV muons and 1.8 x 10 21 useful decays of 20 GeV antimuons, with 
a 50 kt detector at L = 1800 (top), 2900 (center), and 4000 km (bottom), with bm\ x = 5 x 10~ 5 
(left) and 10 -4 eV 2 (right). The solid and long-dashed curves correspond to the CP conserving 
cases 5 = 0° and 180°, and the short-dashed and dotted curves correspond to two other cases that 
give the largest deviation from the CP conserving curves; along these curves sin 2 29i% varies from 
0.0001 to 0.01, as indicated. 
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FIG. 12. 3ct error ellipses in the [iV( / u + ), iV(/i _ ])-plane, shown for a neutrino factory delivering 
3.6 x 10 21 useful decays of 20 GeV muons and 1.8 x 10 21 useful decays of 20 GeV antimuons, with 
a 50 kt detector at L = 7300 km, Sm^i = 10~ 4 eV 2 , and 5 = 0. Curves are shown for both signs 
of <5tb 2 2 ; sin 2 26^ varies along the curves from 0.0001 to 0.01, as indicated. 
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